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Wave function statistics at the symplectic 2D Anderson transition: bulk properties 
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The wave function statistics at the Anderson transition in a two-dimensional disordered elec- 
tron gas with spin-orbit coupling is studied numerically. In addition to highly accurate exponents 
(ao=2. 172 ± 0.002, T2=l. 642 ± 0.004), we report three qualitative results, (i) The anomalous dimen- 
sions are invariant under q — > (1 — q) which is in agreement with a recent analytical prediction and 
supports the universality hypothesis, (ii) The multifractal spectrum is not parabolic and therefore 
differs from behavior suspected, e.g., for (integer) quantum Hall transitions in a fundamental way. 
(iii) The critical fixed point satisfies conformal invariance. 

PACS numbers: 72.15.Rn, 05.45.Df 



Disordered electron systems that are confined to two 
spatial dimensions (2D) cannot support a true metallic 
state because of Anderson localization^ The underlying 
physics relates to an interference-enhanced return prob- 
ability of quantum mechanical particles due to repeated 
backscattering of the same (quenched) disorder configu- 
ration. There are exceptions to the rule, however. For 
instance, if spin-orbit scattering exists, the return prob- 
ability is not enhanced but even depleted and the metal- 
lic state survives^ Universal properties of such metals 
are described by the symplectic symmetry class of Gaus- 
sian random matrix theories. By increasing the disorder 
strength W, a metal-insulator (i.e. Anderson transition) 
can be driven in these materials. Its universal properties 
have been studied intensively in the last two decades. 

One of the controversial questions in the late 1990s 
concerning the symplectic transition in 2D was about the 
numerical value of the critical exponent v that describes 
the divergence of the localization length when the disor- 
der approaches its critical value: £ ~ \W — W c \~ u . In 
recent work, Asada et al. have made a very convincing 
case in favor of v—2.75 (overview in Table [T| employing 
the SU(2) model. 3 A work by Markos and Schweitzer- 
comes to a similar conclusion, v « 2.8 ± 0.04, within the 
Ando model and the debate is now settled. 

However, this latter work not only has helped to fix v, 
it also has reemphasized that another important topic is 
still unresolved. Recall that the critical wave functions, 
^(x) at the boundary between insulator and metal obey 
a multifractal statistics. 5 This implies that the moments 
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scale with system size L, introducing the exponent spec- 
trum T q . (The angular brackets denote a combined spa- 
tial and ensemble average.) A precise numerical determi- 
nation of T q has not been undertaken yet. The numerical 
work presented in this letter is an attempt to close this 
gap. 

There are several good reasons why one would like to 



scrutinize the nature of T q more closely. For one thing, 
the wave function statistics can be measured, in principle, 
and promising steps in this direction were made not long 
ago.— 

But also important questions concerning our con- 
ceptual understanding of the localization-delocalization 
transition are closely related to multifractality. First, the 
analytic structure of r q is a specific characteristic of the 
critical field theory of the transition describing scaling of 
the local density of states. For example, it has been pro- 
posed that the (integer) quantum Hall transition exhibits 
reduced anomalous dimensions S„, 



d{q-l) + 5 q q(l-q), 



(2) 



with a special property: 5 q does not depend on q, such 
that T q is parabolic and also invariant under q — » 1— q. 
Very recently, it has been predicted 7 - - based on exact 
results for the nonlinear a model and invoking the uni- 
versality hypothesis - that this last symmetry is a general 
property of all transitions belonging to the conventional 
Wigner-Dyson classes. That is 
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(3) 



should hold. A numerical verification beyond the frame- 
work of the power law random banded matrix model has 
not been reported yet. This would be an interesting test 
of universality, since it does not only rely on compar- 
ing quantitative values for some few exponents - which 
has been the usual procedure - but rather refers to the 
analytic structure of an exponent spectrum. Note that 
Eq. ([3]) does not generally hold outside the conventional 
symmetry classes. The spin quantum Hall effect is an 
example for a transition in a nonstandard universality 
class, where Eq. © is manifestly violated^ 

Second, lately it has become clear that near bound- 
aries multifractality differs from the bulk: flat interfaces 
support their own "surface" spectrum r|; in the presence 
of corners yet another spectrum is superimposed, etcJ^ 
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5i=0.16 ± 0.02 
5 2 =0.185 ±0.01 
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1.83 ±0.03 






2.51 ±0.18 
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TABLE I: Overview of results for the symplectic transition in 
model 2 ®; MAt (MAa): multifractal analysis based on scaling 
SU(2) model 1 ®; S q : reduced anomalous dimension, see Eq. ([2]) 
with the latest work. 

Also, in principle, an edge could break a bulk symmetry 
and thus would not even share the bulk universality class. 
In fact, the unraveling of surface multifractality could 
lead to a paradigmatic shift of our present understanding 
of critical wave function statistics. Clearly, a prerequisite 
for all this is a detailed knowledge of bulk properties. 

Third, finally, a relation between Sq and the ratio A c 
of width and localization length of quasi- ID strips exists: 

A c = 1/tt5 , (4) 

which is exact if the critical 2D fixed point is conformally 
invariant^ It is believed that conformal invariance (CI) 
is a generic property of localization-delocalization tran- 
sitions in 2D. For instance, it has been demonstrated to 
hold at the integer quantum Hall transitiom 12 i 13 Excep- 
tions are not known so far, but Eq. ((4]) can be used 
as a test of CI. In this respect, recent numerical results 
are alarming. It is reported^ that <5o=0.107 ± 0.005 and 
A c =1.87 ± 0.02; thus the product ttA c 5 = 0.629 ± 0.036 
would signal a strong violation of Eq. Q and therefore 
absence of ClJ^ 

In this Rapid Communication, we present a numerical 
high-precision study of 5 q at the 2D-symplectic transi- 
tion. Our particular aim is to answer three qualitative 
questions, (i) Is 5 q a constant, so r q is parabolic? (ii) If 
not, does it obey the symmetry relation Eq. ([3]) confirm- 
ing the universality hypothesis? (iii) Is the fixed point 
conformally invariant? 

Most earlier works analyzed typical moments in small 
ensembles, where finite-size effects make it difficult to 
obtain reliable error bars. By contrast, we employ scaling 
of typical and average moments in very large ensembles 
with big system sizes. Errors can thus be reduced by 
almost an order of magnitude. In order to cross-check, 
we analyze the two most important microscopic models. 
Results thus obtained agree very well. Specifically, we 
find that 6 q is not a constant and the symmetry relation 
d3j) is satisfied. 

On a quantitative level, we obtain ^2=0.180±0.002 
(both models), <5 =0.173±0.003 (Ando model), and 
(5 =0.172±0.002 (SU(2) model). Together with Eq. @D 



two dimensions. AM: Ando model-*; EZM: Evangelou-Ziman- 
of typical (average) amplitudes; TM: transfer matrix; SU(2): 
Entries for the same model are in chronological order, starting 

and the earlier result 3 A c = 1.843 we arrive at ttA c Sq = 
0.996±0.012. Thus numerical evidence is provided that 
the symplectic fixed point obeys CI, in agreement with 
general expectations. 

Models: We consider a tight-binding Hamiltonian on 
a two dimensional square lattice with nearest neighbor 
coupling 

where c] (c { ) denotes a creation (annihilation) opera- 
tor of an electron with spin a on site i. 

In the Ando model^, the on-site energies ej are taken 
independently from the interval [— W/2, W/2] with a ho- 
mogenous distribution. The hopping matrix Vi ttT -j, a > re- 
flecting the spin-orbit coupling is chosen as 

Vi <a -i + k, a > = (Vb exp(i6 k o-k)) a ,cT> > k = x,y, (6) 

with o~ x ,o~y denoting Pauli matrices and the parameters 
Vq = 1 and 9k = n/6. We have determined the critical 
disorder strength independently via analysis of the crit- 
ical level statistisJ^ Our finding W c =5.85±0.025 agrees 
well with earlier worki^ 

The second model, the SU(2) model, has been intro- 
duced by Asada, Slevin, and Ohtsuk;! 8 . In addition 
to the on-site energies ej, now also the hopping matrix 
Vi ja -j t a-> is random. It is taken to be uniformly distributed 
over the entire group SU(2) using the group invariant 
(Haar) measured 

H is implemented on square L x L-size lattices with 
periodic boundary conditions. For our numerical diag- 
onalization of the resulting 21? x 2L 2 matrices we use 
an inverse iteration routine coupled with direct sparse 
solvers in order to obtain the eigenvalues and wave func- 
tions with energies closest to zeroJ^ (Cf. Ref. |20| .) 

Multifractal analysis: Our multifractal analysis pro- 
ceeds by analyzing the scaling behavior of the average 
moments of wave function amplitudes, Eq. (fTj). 

In order to analyze the critical behavior we take the 
disorder value W c =5.84 (for states at energy zero being 
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FIG. 1: Reduced anomalous dimension 8 q as defined in Eq. 
© for the Ando model (dashed, 14^=5.84) and the SU(2) 
model (solid, W / c =5.953). Additionally, anomalous dimen- 
sions 5 q obtained from typical inverse participation ratios are 
shown (o) for the latter model. Dashed lines indicate the 
estimated error (2a) in So- Inset: blowup of the solid line 
behavior near q=0.5 now represented by o. Data near q—0 
and q=l suffer from noise amplification (dividing by q(l— q) 
in Eq. (J5J) and have therefore been omitted. Filled symbols 
(•) show original trace after reflection at q=0.5. Dot-dashed 
line indicates parabolic fit (offset: 10 -3 ) with 5 1/ / 2 =0.1705 
and curvature <5"/ 2 =0.0043. 



critical) in the Ando model. For the SU(2) model we 
employ W c — 5.953 in order to have a mobility edge at en- 
ergy e=l£ The average ([T]) has been performed over an 
ensemble of wave functions that have been calculated in 
systems of sizes £=16, 24, 32, 48, 64, 96, 128, 192, 256 (the 
last two values were not used in all cases). For each dis- 
order realization 64 wave functions closest to the critical 
energy have been taken into account; all together the 
number of wave functions in the ensemble is typically 
4 x 10 7 (£=16) to 3 x 10 5 (£=256). 

The exponents r q are readily extracted from a power- 
law fit as suggested by Eq. (TT]). 21 In Fig. Q]we plot the 
reduced dimensions S q defined in ^ as obtained for both 
models. It incorporates our three main results. 

(i) We determine <5 =0.172 ± 0.002. The value satisfies 
Eq. (|2]) and thus the consistency check on CI is positive. 
The good accuracy stems mainly from large statistics and 
the fact, that finite-size corrections in the SU(2) model 
turn out to be extremely small at q < 1.5. As can also be 
seen from Fig. [TJ the Ando model gives a similar result. 

(ii) The function S q satisfies the symmetry relation Eq. 
([3]). Thus the universality postulate is confirmed. The 
inset of Fig. [T] shows that part of the full curve 6 q , for 
which numerical data are available at both points, q and 
its image 1 — q. (The numerical procedure that we work 
with is limited to q> — 1; more negative values would 
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FIG. 2: Test function £l q (L) highlighting variability of S q with 
q for two q values. Solid lines represent power law fits, with 
a slope representing S q — #1/2; for values, see legend. Slight 
deviations between models are due to larger errors in finite- 
size extrapolation of Ando model. For that model, results for 
two values of W are given to illustrate that the uncertainty 
in W c is not a precision- limiting factor. Note that £3/2— <5i/2 
agrees well with curvature 8'{/ 2 seen in Fig. [1] 



require a coarse graining in order to overcome the diver- 
gence of the moments (1) related to zeros of the wave 
functions.) A symmetric shape of the curve is clearly 
displayed in the regime of best accuracy, —0.5 < q < 1.5. 

(iii) The set of exponents S q does not reduce to a con- 
stant, e.g., S q has a small but nonzero curvature S'^ 2 . 
Detecting <5"y 2 requires high-precision data, because the 
numerical window is limited to q < 2.0. At larger values, 
(a) finite-£ effects proliferate (in Ando model faster than 
in SU(2)), so deviations between solid and dashed lines 
increase. And (b) moments ((\^\ 2q )) for large q probe the 
tails of the distribution function, so that typical values 
and averages differ from each other. Then, error bars 
tend to become large due to undersampling3 The part- 
ing of the three curves at q > 2 visible in Fig. [1] is a 
consequence of these effects. 

As a sensitive test for variability of S q we investigate 
in Fig. [5] the ratio 
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encompassing only unprocessed data. It scales as 
n q (L)~L- 5 i +5 i/ 2 and therefore any slope in In f2 signal- 
izes that S q deviates from 5 1 / 2 =0.1705 ± 0.001. Data for 
il q at <7=1.5, 2.0 are shown in Fig. It clearly exhibits 
a linear trace with the nonzero slope indicative of curva- 
ture in Sq. Note that finite-size effects are very small, so 
that S q — S 1 / 2 can be extracted with good accuracy. 

A more conventional object than S q to characterize 
the wave function statistics is the Legendre transformed 
f(a) — qa — r q , a q =dr/dq, displayed in Fig. [3) Even 
though we have obtained r q only for q> — 1/2 and there- 
fore are restricted to a<ai/2, the spectrum can be re- 
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FIG. 3: f(a) spectrum from data of Fig. \T\ SU(2) model. 
f(a) is slightly asymmetric and not a parabolic function, 
which would have meant f(a)—2 — (a — 2 — So) 2 /4<5q. 



constructed also at values a>a\/2 by making use of Eq. 
Q. 7 Then deviations from parabolicity obtrude. 

Summary: The multifractal spectrum of wave func- 
tions at the 2D symplectic Anderson transition has been 
calculated in the Ando and SU(2) models with high pre- 
cision. On a qualitative level, our results demonstrate 
that the critical fixed point is conformally invariant with 
a nonparabolic spectrum r q . Furthermore, 5 q =5i- q , as 
predicted from calculations within the nonlinear a model 
and thus supports the universality hypothesis. 

We thank L. Schweitzer and K. Yakubo for useful cor- 
respondence and A. D. Mirlin for valuable discussions 
and suggestions on the manuscript. While finalizing the 
manuscript, we learned about a closely related project, 
with partly overlapping results.— 
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